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Abstract 

Cross sections for hyperon pair production from two-photon collisions, 77 — > AA, S^S", are 
measured in the 2 — 4 GeV energy region at Belle, using 464 fb~^ of data. A contribution from the 
intermediate resonance ?7c(l'S') is observed, and the products of the two-photon width of the rfc{lS) 
and its branching ratios to AA and S^S^ are measured. The results will help test QCD models. 

PACS numbers: 12.38.Qk; 13.60.Rj; 14.40.Gx 
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INTRODUCTION 



Measurement of hyperon-pair production from two-photon colhsions is important for 
the study of QCD models and flavour symmetry. Previous measurements 0, 01 have been 
performed with very limited statistics; the results of different experiments are not consistent. 
General theories of hard exclusive processes in QCD 0, 0, and the pure-quark picture 
result in an expected cross section for baryon-pair production from two-photons that is 
one-order of magnitude below the experimental data. In order to explain the experimental 
observation, various models were proposed. For example, the diquark model 0, H, |^, 
achieves better agreement in the absolute size of the cross section. In handbag approaches 



the process is factorized into a hard 77 — >• qq sub-process and a soft qq — >• baryon- 
pair transition; the expected cross section is determined by the so-called annihilation form- 
factors. In these calculations, all effects from flavour symmetry breaking are neglected. More 
recently, predictions were also provided within pole and resonance approaches 3 • 



Using 464 fb~^ of data collected at Belle, hundreds of events have been obtained for each 
of the two channels 77 AA, E'^E°. Contributions from r]c resonances are also observed. 
The measured cross sections will help to test all the existing models more precisely. 



BELLE DETECTOR AND EVENT SELECTION 



Data are recorded with the Belle detector |13| at KEKB 
e~^e~ collider operated at 10.58 GeV center-of-mass (cm.) energy. The following Belle sub- 
systems are of importance for our analyses: the central drift chamber (CDC), the aerogel 
Cherenkov counters (ACC), the time-of-flight scintillation counters (TOF) and the CsI(TZ) 
electromagnetic calorimeter (ECL). The CDC measures the momenta of charged particles 
and provides precise (6%) dE/dx measurements. The ACC measures the number of pho- 
toelectrons produced by highly relativistic particles. The TOF measures the time of flight 
of particles with a 100 ps timing resolution. The ECL detects photons and the deposited 
energy of particles with a resolution oi ge/ E = 1.5% (2.0%) at 1 GeV (0.1 GeV). 

In processes with quasi-real two-photon collisions e^e" —>■ e~*"e~7*7* — > e"*'e~ AA(E°E°), 
the scattered electrons go down the beam pipe and thus only the produced hyperon-pair 
can be detected. The 77 axis can be approximated by the beam direction in the e~^e~ cm. 
frame. Taking into account the fact that A(A) can be reconstructed from its pTT~{7r^p) 
decay products and S°(S°) from A7(A7), candidate events including pn^n^p candidates 
are searched for in a low-multiplicity data stream, where exactly four tracks, two positively 
and two negatively charged, are required. The sum of the magnitudes of the momenta 
of all tracks and the total ECL energy are restricted to be below 6 GeV/c and 6 GeV, 
respectively. Each of the tracks satisfies the first level conditions: pt > 0.1 GeV/c, dr < 5 
cm and \dz\ < 5 cm, or two or more tracks satisfy the second level conditions: pt > 0.3 
GeV/c, dr < 1 cm, \dz\ < 5 cm and 17° < 6 < 150°. Here pt is the transverse momentum 
with respect to the beam axis, dr and dz are the radial and axial coordinates of the point 
of closest approach to the nominal collision point respectively, and 6 is the polar angle of 
momentum with respect to the electron beam. Among the first level tracks, two tracks 
with opposite charge are required; using the second level tracks, the invariant mass and the 
missing mass squared (with a zero mass assumption) have to be smaller than 4.5 GeV/c^ 
and larger than 2 GeV^/c^, respectively. 

Backgrounds are further reduced by a particle identification (PID) algorithm; two tracks 



14 1. which is an asymmetric 
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FIG. 1: Two-dimensional distribution of the normalized likelihood of the p and the p candidates, 
for the events passing all selection criteria up to |SpJ'(AA)| < 0.2 GeV/c except the cuts on the 
normalized likelihood. In (a), the dark bins show the events passing the cuts on the normalized 
likelihood; in (b), the dark bins show the events with either of the two tracks satisfying Lp/L^ < 0.5, 
where x is K, vr, n or e. 

with opposite charge have to pass the p{p) identification, and the other two have to pass the 
vr^ identification. The p{p) identification includes the following conditions : 

1. the difference between the measured and the expected CDC dE/dx is less than 6 times 
the resolution: 

-A{dE/dxY^ 



XaE/cLx 



dE/dx 



< 6^; (1) 



2. the ratio of the associated ECL energy to the momentum is less than 0.9 for the 
positively charged track; 

3. the number of the photoelectrons in the ACC associated with the track is less than 
10; 

4. the likelihoods for each particle assigment are combined to determine the normalized 
likehhood, 

Xp = ^ , (2) 

Lp + Lk + Lt^ + + Le 

which has to be larger than 0.95 (0.2) in case an associated TOF hit is available (not 
available). Each likelihood 



exp 



2(^dE/dx + Xt) 



(3) 



is calculated using CDC {dE/dx) and TOF (time of flight T) information [15i] . 
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FIG. 2: Distributions of the invariant mass of (a) p,iT~ and (b) tt^,P for the events passing all 
selection criteria up to |Epj(AA)| < 0.2 GeV/c except the cut |Amy\^^^)| < 10 MeV/c^, indicated 
by the arrows. The Monte Carlo is scaled by a factor corresponding to the normalization to data 
in the p,7T~ invariant mass range 1114.8 — 1116.3 MeV/c^. The two-dimensional distribution is 
shown in (c). 

Figure 1 shows the two-dimensional distribution of the normalized likelihood for the p,p 
candidates before the cuts on Ap, where the dark region in (a) shows the events passing the 
cuts, and that in (b) shows the events with at least one of the tracks satisfying Lp/L^ < 0.5 
with X being K, ir, /z or e, which is referred to as the backgound-dominant region and is 
almost rejected. 
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For n"^ identification, only dE/dx is used; the pion candidates must have dE/dx consistent 
with the pion hypothesis {\dE/dx — {dE / dx).„±\ < 6 adE/dx) and satisfy a proton veto 
{\dE/dx - (dE / dx)p(^p)\ > 6 adE/dx)- 

The tracks passing the above PID conditions are paired into pn^i^n^p) combinations, 
and each of them has to pass the A (A) vertex reconstruction. The difference between the 
invariant mass of pTT^^TT^p) and the A(A) mass, Am^j^-j, is required to be less than 10 
MeV/c2 (Fig. 2). 
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FIG. 3: Distribution of AA invariant mass for (a) the events passing the selection criteria up to 
|$]pj(AA)| < 0.2 GeV/c in the 2.5 — 3.4 GeV/c^ mass region, and for (b) the events passing the 
exclusive 77 AA selection with |Spj(AA)| < 0.05 GeV/c. In (a), the data in the mass range 
3.08 — 3.11 GeV/c^, corresponding to the region of radiative return to J/jp background, are not 
used in the fit. 



In order to reject most of the non-exclusive 77 - 
verse momentum balance in the cm. frame of e~^e~ 



* AAX, S'^S^X backgrounds, the trans- 
beams is required to satisfy 



isp:(AA)i 



EC 



< 0.2 GeV/c, 



(4) 



where p^* denotes the transverse momentum of each track in the e~^e cm. frame. Disregard- 
ing the two photons from the S°S*^ A7A7 decay, the |Ep^(AA)| peak of the 77 E°E° 
events is around 0.1 GeV/c according to Monte Carlo simulation, so that most of them 
pass the condition in Eq. (4). In total, 1628 events survive after all of the selection crite- 
ria above. At this stage, peaks around 2.98 and 2.80 GeV/c^ due to the intermediate rjc 
resonance, corresponding to the process "yy ^ f]c ^ AA, S'^S°, are observed in the distribu- 
tion of AA invariant mass (Fig. 3(a)). Details will be given in the next section. Exclusive 
77 — > AA events are further required to satisfy |Sp^(AA)| < 0.05 GeV/c, where 549 events 
are obtained (Fig. 3(b)). 
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FIG. 4: Distributions of the invariant mass of (a) A, 7 and (b) A, 7 for the events passing all 
selection criteria for exclusive 77 — > Tj^T,^ except the cut |Am20(20)l < MeV/c^ indicated by 
the arrows. The Monte Carlo includes the backgrounds from 77 — > AA, AAtt'^, T,^T,^tt^ and is scaled 
by a factor corresponding to the normalization to data in the A, 7 invariant mass range 1181 — 1201 
MeV/c^. The two-dimensional distribution is shown in (c). 

In the sample of 1628 events satisfying |Sp*(AA)| < 0.2 GeV/c, exclusive 77 — >■ S°E° 
events are selected as follows. Each A (A) is paired with a photon candidate, where photons 
are selected from the ECL clusters satisfying the following conditions: 

1. the cluster is not associated to any CDC track; 
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2. the total energy of the cluster is between 50 and 200 MeV; 



3. the cluster is photon-like and isolated: Eg/ > 0.9, where Eg and E25 are the energy 
deposition in the 3x3 and 5x5 matrix around the ECL crystal with the maximum 
energy. 



< 



The A7(A7) pair has to satisfy lAm^o^j^o-) 

nations the one with the smallest |Amj.o(so)| is selected (Fig. 4). Here |Am20(so) 
difference between the invariant mass of A7(A7) and the S°(S° 



100 MeV/c^, and among all possible combi- 

is the 



mass. 



o 



O 

(N 



CZ3 



> 



5-1 



18 
16 
14 
12 
10 
8 
6 
4 
2 




ll 



1 



I 



Fit 




2.6 2.8 



3.2 3.4 3.6 3.8 



invariant mass (GeV/c^) 

FIG. 5: Distribution of S^S" invariant mass for the events passing the exclusive 77 
selection requirements. 



In order to remove the background from exclusive 77 AA events, events with 
|Sp*(AA)| < 0.06 GeV/c2 and |Sp*(A7A7)| > 0.02 GeV/c^ are excluded, where |Sp*(A7A7)| 
is the transverse momentum balance including the contribution from the two photons. Fur- 
thermore, in order to reject other non-exclusive backgrounds, |Ep^(A7A7)| and |Sp^(AA)| 
are restricted to below 0.14 and 0.16 GeV/c^, respectively. About 11% of the events have 
two or more photon candidates, but the A and A are not paired with two different pho- 
tons; these events are simply omitted. After all the above requirements, 212 events remain 
(Fig. 5). 



OBSERVATION OF 77 ^ 7?c ^ AA, SOSO 

From the sample of the 1628 events satisfying |SpJ'(AA)| < 0.2 GeV/c, peaks around 2.98 
and 2.80 GeV/c^ in the distribution of the AA invariant mass (Fig. 3(a)) are identified as 
decays of the rjc{lS) resonance [l^. The former corresponds to the •j'j rjc ^ AA signal; 
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while based on a Monte Carlo study (Fig. 6), the latter is identified as the effective signal 
from 77 — r/c — ^ E°E° — A7A7, where the two photons in the final state are not detected. 
When the exclusive E'^S" selection requirements are applied, the latter signal is recovered 
at the mominal t]c{1S) mass in the A7A7 invariant mass distribution (Fig. 5). 
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FIG. 6: Distribution of AA(S°i;0) invariant mass from a Monte Carlo 77 
passing the AA(S''S'') selection, respectively. 



S^SO sample 



The rjc — >■ AA signal in Fig. 3(a) is fitted with a Breit-Wigner function, and an rjc yield of 
101.2 ± IGAtll events is obtained. The Breit-Wi gner function for the rjc signal is smeared 
with detector resolution determined by Monte Carlo simulation, and the mass and the width 
of the ?7c are also free parameters in the fit. The t]c mass and width obtained from the fit are 
consistent with those reported in Ref. Backgrounds are fitted with a smooth function 
(exponential of a first-order polynomial) for the continuum part and a Gaussian for the 
E°EO contribution. For the latter, the Gaussian mean value and the width are also 

> AA 

where Ls 



Vc 

free parameters in the fit. We obtain a statistical significance of 6.60" for the r], 
signal from the fit. Here the statistical significance is defined as 



2 HLo/Ls 

and Lq denote the maximum likelihoods of the fits with and without a signal component, 
respectively. A similar fit is performed for the rjc — > S'^SO signal in Fig. 5, using the same 
type of functions as above for the signal and the continuum background. An rjc signal of 
36.1 ± 9.2~^i'^ events with a statistical significance of 3.9a is obtained. For the systematic 
error in the t]c yield, effects from the uncertainties of the continuum background shape and 
the energy detection resolution of signal are taken into account. 

The rjc yield (A^rjJ is converted to a product of the two-photon width of ric and its 
branching fractions to A A or according to the formula 



r^^(r/,)x5(r/,-.AA,S°S0) 



47r2 D e U^t dL^^/dW. 



(5) 



77 
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using the luminosity function dL^^/dW^y (Eq. (7)) determined at the energy of the rjc mass 
{rrir^J with efficiency e from Monte Carlo. Here D denotes the branching fraction for AA 
p7i~n^p, which is equal to 40.83 ± 0.64% |16)], and Lint denotes the integrated luminosity 
which is equal to 464 fb~^. The results are shown in Table I, where the systematic errors 
are estimated taking into account the sources listed in Table II. The value of B{ric — >■ AA) 
given in Table I is in agreement with the Belle measurement from B meson decays [l7| . 



TABLE I: Measured T^^{ric) x B{r]c AA,S°SO) and B{r]c AA,E°SO). The latter is obtained 
using r^^(77c) = 7.0jiQ'g keV Q|. The errors are statistical and systematic respectively. The last 
error in B[r]c — > AA, S'^S'^) is from the uncertainty in T^^{r]c)- 



rjc decay mode 


r^^(r/e) X B{r]c ^ AA,SOSO) (eV) 


B{r]c AA,i;°SO) X 10^ 


rjc AA 
r]c SOSO 


6.21 ± 1.01 toil 
9.80 ± 2.50 till 


0.89 ± 0.14 

1.40 ± 0.36 t'o\t 



TABLE II: Systematic errors (%) for the measured T^^{r]c) x B{rfc KK^TPTP) and B{ric — > 
AA, E^EO) 



Source 


Vc AA 


^ S^EO 


Integrated luminosity 


1.4 


1.4 


Luminosity function 


4 


4 


Branching fraction of AA pTT~'ir~^p 


1.6 


1.6 


Monte Carlo statistics 


1.0 


2.4 


Trigger efficiency 


4.3 


4.3 


Particle identification efficiency 


4.5 


4.5 


Photon selection efficiency 





5.9 


Background shape and energy detection resolution of signal 


+1.2 
-3.0 


+0.0 
-3.3 


Total 


+7.9 
-8.3 


+10.0 
-10.5 



MEASUREMENT OF THE CROSS SECTIONS FOR 77 ^ AA, E°EO 

The cross sections for 77 — > AA(E°S°) are measured, using the data sample of 549 (212) 
events (Fig. 3(b) and 5) passing the exclusive 77 — >■ AA(S°E'^) selection requirements. The 
number of events AN(W^-y, \ cos6'*|) and the efficiency e(W^^, \ cos6'*|) are determined for 
two-dimensional bins of W^^ and | cos^*|, where W^^ and 6* are the two-photon cm. energy 
and the cm. angle of the hyperon, respectively. The background from radiative return to 
J/ip contributes most of the excess events in the range of W^^ = 3.08 — 3.10 GeV, and 
thus the data in that narrow bin are not used. The ratio AN/e is then converted to the 
differential cross section, according to the formula 

ti^^^^AA,SOSo(^77) _ AA^/g 

d\cose*\ D Lint §9^ AW^^ A|cos^*| 
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The luminosity function ^r^, as a function of W^^, is defined by 

^77^AA(S0S0)l'^77J dW^^ [7) 

and is calculated by TREPS [l^ using the equivalent photon approximation method jl8l |. 
The effects from longitudinal photons are neglected. For event generation in TREPS, the 
maximum virtuality of each of the two photons, Qf and Ql, is limited to 1 GeV^, and a form 
factor term is introduced for the high-Q^ suppression effect, {1 + Q'l/W?^^^)~'^{l + Ql/W^^)~'^. 
The systematic uncertainty in the luminosity function is estimated by comparing TREPS 
to a QED calculation including all order diagrams j20|, and an agreement within 3 — 5% 
was reported for W^^ = 2 - 4 GeV pj, lU 



Overall detection efficiencies from Monte Carlo {e) are obtained from the signal samples 
generated by the TREPS codes Other samples of the type 77 AAvr^, S'^E'^vr^ are 

generated using the GGLU [2^ code in order to study non-exclusive backgrounds. Detector 
simulation is based on GEANT3 Typical values of e from Monte Carlo simulation in 
the region with | cos 61*1 < 0.1 and W^^ = 2.5 - 4.0 GeV are ~ 2% - 9% for 77 AA and 
~ 0.3% - 3.5% for 77 S^SO (Fig. 7). 




FIG. 7: Overall detection efficiency as a function of W^^ and | cos^*| for (a) 77 AA and (b) 
77 — > S'^S'^ events. 

Knowing the differential cross section, the total cross section can be obtained by summing 
over I cos 6* I : 



cr. 



77-AA,SOEo(W^77) = L d\cOs9*^ A| COS0 | (8) 



which is performed with A|cos^^*| = 0.1 for |cos6'*| up to 0.6. The results are shown 
in Fig. 8, where the contribution from the intermediate rj^ decay is included. Previous 
measurements from CLEO ^J, L3 and the 77 —>■ pp cross sections from Belle jl^, as 
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FIG. 8: Measured cross sections for 77 AA,T,^T,^,pp. The error bars are purely statistical for 
Belle, and the contribution from the intermediate rjc decay is included. The theoretical curves 
shown are the predictions from the diquark model j9|. 



well as the predictions given by the diquark model are also shown for comparison. The 
major sources of systematic error are listed in Table IIL 

The Monte Carlo trigger efficiency for the low energy trigger is checked using experimental 
data. Over 80% of the signal events pass this trigger. The low energy trigger requires a 
0.5 GeV ECL total energy sum and at least two tracks that go through three inner CDC 
trigger-layers. The trigger efficiencies of the two sub-parts of the low energy trigger (i.e. 
the ECL and the CDC parts) are separately obtained using the experimental data passing 
the two-track trigger and the high energy trigger respectively, where the former requires 
at least two tracks that go through all CDC layers and the latter is based on a 1 GeV 
threshold for an ECL total energy sum 3, 2^- The overall trigger efficiency is ~ 60 — 80% 
corresponding to an average transverse momentum of the A and A from 0.3 to LO GeV/c 
data. Corrections for the Monte Carlo trigger efficiency are implemented according to the 
data, and the systematic error is 4.3%. 

The systematic uncertainty of the event selection efficiency is predominantly coming from 
PID and the photon selection. The latter denotes the S^E^ selection up to the paring of 
A7A7 to which is applied to the events with an identified AA pair. The overall PID 

and the photon selection efficiencies are ~ 90% (~ 86%) and ~ 55% (~ 45%) at W-y^ = 3 
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TABLE IIL Systematic error (%) for the measured cross sections for 77 AA, S'^S*' 



Source 


77 aa 


77 ^ SOso 


Integrated luminosity 


1.4 


1.4 


Luminosity function 


5 


5 


Branching fraction of AA pir^n^p 


1.6 


1.6 


Monte Carlo statistics 


0.5 - 8.6 


1.0 - 8.1 


Trigger efficiency 


4.3 


4.3 


Particle identification efficiency 


4.5 


4.5 


Photon selection efficiency 





5.9 


Residual background 


15 - 26 


10-37 


Total 


17-29 


14-39 



GeV (4 GeV). The accuracy of the Monte Carlo efficiencies are checked from the rjc yield 
in the data, where a realistic PID or photon selection efficiency can be estimated from 
the number of events with and without the application of the PID or the photon selection, 
respectively. The systematic errors are 4.5% and 5.9% corresponding to PID and the photon 
selection efficiency respectively, for W^^ in the 2.7 — 3.2 GeV region. 

One of the major backgrounds is non-exclusive decays surviving the transverse momentum 
balance cut, and the contamination from such backgrounds is checked from the Monte Carlo 
samples of 77 AAvr^, E°S°7r°. The most conservative estimates are made assuming these 
backgrounds are the same size as the signals. Another significant residual background type 
comes from the signals themselves, 77 — >■ AA, S°S°, which could contaminate each other 
according to a Monte Carlo study. A possible background from 77 — » AS°(S°A) is omitted 
based on the expected level from current theoretical predictions j9|,|ll|. The backgound from 
radiative return is at the level of 1% or less. The overall residual background is 15 — 26% 
for 77 —>■ AA and is 10 — 37% for 77 — > E^S^^, where the ranges are determined by using 
the least and the most conservative estimates for the non-exclusive background in the data. 
Since it is difficult to subtract all these possible backgrounds accurately, the estimated values 
above are assigned as systematic uncertainties in the measured cross sections. 

CONCLUSION 

Using 464 fb~^ of data, cross sections for hyperon pair production 77 — > AA, S'^SO have 
been measured at Belle. Intermediate rj^ resonances are observed, and the products of the 
two-photon width of the 7]^ and its branching ratios to AA, are determined. The result 
for rjc S^S'^ is the first measurement. 

The measured cross sections for 77 AA from Belle are compatible with L3 P] results. 
The data show that the cross sections for the three channels 77 AA, E°E°,j9p converge 
at high energy. The deviation of the theoretical predictions from data implies that current 
models have to be improved; the effects from flavour symmetry breaking could be non- 
negligible in the measured energy range. 
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